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Abstract 



The early universe is believed to have undergone a QCD phase transition to hadrons at about lOps 
after the big bang. We study such a transition in the context of the non-detailed balance Horava-Lifshitz 
theory by investigating the effects of the dynamical coupling constant A in a fiat universe. The evolution 
of the relevant physical quantities, namely the energy density p, temperature T, scale factor a and the 
Hubble parameter H is investigated before, during and after the phase transition, assumed to be of first 
order. Also, in view of the recent lattice QCD simulations data, we study a cross-over phase transition of 
the early universe whose results are based on two different sets of lattice data. 

o 
cr 

vL| ■ 1 Introduction 

bJQ' 

Recently, a power-counting, ultra-violet {UV) renormalizable quantum gravity theory has been proposed [TJ__], 
' stirring a considerable amount of interest in its wake. This theory is principally based on a scalar field theory 

, introduced by Lifshitz [3] to explain quantum critical phenomenon in condensed matter physics. In this theory 

' the action is scale-invariant under the transformations 

VD 

"3". t^£ z t, x l ^lx\ (1) 

in . 

with z = 2, where z represents the dynamical critical exponent. However, the theory proposed by Horava 
has the scaling dimension z = 3 and is often referred as Horava-Lifshitz (HL) gravity and is fundamentally 
non-relativistic. This theory does not have the full diffeomorphism invariance of general relativity, admitting 
only a local Galilean invariance. Einstein general relativity may then be considered as emerging as an infra-red 
(IR) fixed point. Invariably, such a theory may have interesting ramifications when applied to cosmology and, 
not surprisingly, a great amount of interest has been generated in this regard. In particular, study of the early 
universe in the context of the HL gravity H] would represent an interesting area of research. A number 
of works have appeared over the past few years studying the early universe in this context. In [4] it was 
suggested that the divergence of the speed of light in the UV region may resolve the horizon problem. The 
phenomenological implications of the HL theory were addressed in [6] . Black hole solutions and related issues 
were studied very actively in [3 [8j __] [10] , some of the theoretical issues and generalizations of the HL model 
were discussed in [TTJ [T^J [T__> Ell [H] ■ Needless to say, a myriad of other interesting topics and issues can be 
found in the literature. In view of the above, the present work is devoted to the study of QCD phase transition 
in the context of HL gravity in the early universe. 

According to standard cosmology, as the early universe expanded and cooled, it underwent a series of 
symmetry-breaking phase transitions where topological defects may have formed. During its early evolution, 
the universe is believed to have undergone at least two phase transitions. The electroweak theory predicts that 
at about lOOGey there was a transition from a symmetric high temperature phase with massless gauge bosons 
to the Higgs phase, where the SU(2) x U{1) gauge symmetry is spontaneously broken and all the masses are 
generated. The second transition, predicted by quantum chromodynamics (QCD), is believed to have occurred 
at a temperature of about 200 Af eV where there was a transition from a quark-gluon plasma to a plasma of light 
hadrons. At about the same energy we expect that the global chiral symmetry of QCD with massless fermions 
is spontaneously broken by the formation of a quark pair condensates. The nature of the QCD transition 
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affects our understanding of the evolution of the universe, see [TB] for further insight. In a strong first-order 
phase transition the quark-gluon plasma supercools before bubbles of hadron gas are formed. These bubbles 
grow, collide and merge, during which gravitational waves could be produced [T7]. Baryon-enriched nuggets 
could remain between the bubbles, contributing to dark matter. The hadronic phase is the initial condition 
for nucleosynthesis, so inhomogeneities in this phase could have a strong effect on its development [T8] . 

The nature of the phase transition involved in the early stages of the evolution of the universe has been a 
source of some debate. A first order phase transition was the basic assumption under which calculations had 
been done in the past. However, detailed lattice QCD calculations performed over the years has now upset the 
balance in favor of a smooth crossover phase transition. It is therefore the purpose of the present work to study 
the formation of hadrons in the early stages of the evolution of the universe based on the assumption that a 
smooth crossover phase transition has been responsible. For the sake of clarity and comparison, we have also 
performed our calculations based on the assumption of a first order phase transition which we will discuss in 
due course. 

Fundamental to the study of any phase transition in the early universe is a knowledge of the equation of 
state (EoS) which plays an important role in the hot, strongly interacting matter [THl [2D] and hydrodynamic 
description of heavy ion collisions [21,. There are a number of models where the equation of state in a first order 
phase transition is used for describing the collective flow in heavy ion collisions. However, as was generally 
mentioned above, lattice QCD suggests that the transition to the de-confined phase is only a crossover one [22] . 
The effects of the EoS on all aspects of the behavior of the system under study is therefore of great importance. 
For example, it is not obvious to what extent the collective flow is sensitive to details of the equation of state. 
It is now well known that in an ideal fluid dynamics scenario, the anisotropy of the proton flow is most sensitive 
to the QCD equation of state [23] . 

Many attempts have been made in the last two decades to calculate the EoS on a lattice, see [24j[25]. Lattice 
calculations are usually done at finite temperatures and at fixed temporal extent N T where the temperature is 
changed by varying the lattice spacing a, T = l/N T a . The thermodynamic observables are calculated through 
the calculation of the trace of the energy momentum tensor (e — 3p) , known as trace anomaly, which can be 
done in the high temperature region rather accurately which is not the case in the low temperature region. 
Some authors have studied an analysis of the trace anomaly within the Hadron Resonance Gas model (HRG) 
[26] . Such studies show that to reproduce the lattice results for the asqtad and p4 actions of the hotQCD 
collaboration, it is mandatory to distort the resonance spectrum away from the physical one, necessary for 
taking into account the larger quark masses used in these lattice calculations. So, to construct a realistic 
equation of state one could use the lattice data. In a number of works lattice calculations show a reasonably 
good agreement with hadron resonance gas (HRG) model predictions at low temperature [27j [28] [29] and in 
other works such distortions are not needed to describe the data [42] . 

The starting point when studying the evolution of the early universe, or in general any cosmological setting, 
are the Friedmann equations. These equations in HL gravity differ from those of the standard 4D cosmology 
in the presence of the parameter A and of the dark radiation term. We expect this deviation from the standard 
4D cosmology to have noticeable effects on the cosmological phase transitions. The first-order phase transition 
have been studied in HL gravity |30| . in the context of brane- world scenarios |31j and in a Randall-Sundrum 
(RS) brane- world model [32 . It would therefore be interesting to study crossover phase transition by using the 
EoS resulting from lattice studies. In what follows, we shall concentrate on the occurrence of phase transition in 
the early universe in the context of HL gravity, using the type of phase transition and EoS suggested by lattice 
studies, that is, assuming a smooth crossover quark-hadron phase transition. We also present a comparison 
between these results and those when the phase transition is considered to be of first order. 

2 Horava-Lifshitz gravity 

Horava theory is constructed on the basic assumption of anisotropic scaling between space and time, i.e. 



In a theory with anisotropic scaling, time and space are fundamentally distinct. It is then natural to use the 
ADM formalism for which the metric is given by 



where gij, N and iVj are the spatial metric, lapse and shift functions, respectively. These variables are dynamical 
variables which under the above scaling, scale as 



t -> ih, 



x l -> Ix. 



(2) 




(3) 



(4) 
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The action of HL theory is decomposed into a kinetic, a potential and the matter part as 

Shl = Sk + S v + S m . (5) 
The kinetic action in the Horava-Lifshitz gravity (HL) is naturally given as 

S k = J dtdx 3 N^^ (l\,,l\' J - XK 2 ) , (6) 

where Kij [K = K l j) is the extrinsic curvature of spatial slices, defined by 

Kij = -^jVt9ii-ViN i -V j N i ), (7) 

with V being the covariant derivative on the spatial slice whose metric gij is used to raise and lower indices. 
The most general potential for z = 3 can be written as 

S v = ^ J dtd a xi/gN[a + iR + -yxR 2 + -/ 2 RijR ij + ^e ijk R i iV j R{ + ^R 3 + a 2 RR ij Rij (8) 

+ a 3 R{R^Ri + <j A V t RV l R + t T B V i R ik V i R ik ]. 

On the other hand, the above ten coefficients can be reduced to three, denoted by /i, Aw, C> using the detailed 
balance condition [2j Q] 

teW gVAw kV(i-4a) _ ^ B y (Q) 

° ~ 16(3A-1) ' ' — 16(3A-1) ' '! — 64(3A-1) ' ' 2 — 16 ' S — 4C 2 ' W 

where k 2 [= and Aw are the Einstein coupling and cosmological constants, respectively, and finally 

S m = J dtdx 3 N^L mi (10) 

where L m — L m (N,Ni,gij,(f>) is the lagrangian density of matter fields, denoted collectively by <f>. In the IR 
limit the HL action simplifies to 

S Sj ~ J dtdx 3 N^f§[a{K i:j K 13 - XK 2 ) +-/R + a}. (11) 

Defining x° — ct and 

A = l, c= x fl, 16ttG= x [^, A eh =~, (12) 
V ct V a 2a 

will reduce this action to that of the Einstein-Hilbert 

Seh = ( dx 4 y/§(R4 - 2 A E h) (13) 



i— / dtdx 3 ^/gN[K l:j K 11 -K 2 + (R- 2 A)}. 
ttG J 



WirG 

The full space-time metric g^ is given by 

3oo = -N 2 + g lJ N % N 3 , g 0l = N i} 

9ij=9ij, V^= N V9- ( 14 ) 
There is a condition inferred from the study of some cosmological models [33] , given by 

a(3A-l)>0. (15) 
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3 Field equations in Horava-Lifshitz gravity 

Now, in order to focus on a cosmological framework, we use the FRW metric 

N=l, 9ij = a 2 (t) 7ij , N* = 0, 

with 

dr 2 



7ij 



dx l dx J = 



1 - kr 



r 2 dV&. 



The homogeneous and isotropic metric is then written as 

dr 2 



ds 2 = -dt 2 + a(t) 2 



1 — kr 2 



+ r\d6 z + sin 2 6 d<j) z ) 



(16) 



(17) 



(18) 



where k — — 1, 0, 1 represents an open, flat or closed universe, respectively. For the matter energy- momentum 
tensor we restrict our analysis to the case of a perfect fluid given in comoving coordinates by 



(19) 



where p, p and are the energy density, isotropic pressure and four-velocity of the cosmological fluid. We 
obtain the gravitational field equations of the HL theory, corresponding to the line element (TTB1) and p7|) as 



3(3A- l)H 2 = — p + 6 



°_ + kj + 2fc 2 (3 7 i + 72 ) 4fc 3 (9gi+3cr 2 + g 3 ) 
6 a 2 a 4 a 6 



(20) 



with 



(3\-1)(H+Ih 2 ) = ^ P - 



a 
6 



3^2 



p + 3H(p + p) = 0. 



2fc 3 (9cri + 3cr 2 + o- 3 ) 



(21) 



(22) 



Since the early universe is known to be flat to within a good approximation, in the follow we shall focus on the 
non detailed balance condition with k = 0. 



4 First order quark-hadron phase transition 

In this section, we present a discussion on the derivation of thermodynamical quantities of the quark-hadron 
phase transition which corresponds to the first order phase transition. It is now generally accepted that such 
phase transition has occurred at the cosmological QCD scale about t ~ 10~ 5 s after the big bang and at a 
temperature T ~ 200M eV where the Hubble radius was about 10 km. Note that the mass inside the Hubble 
volume is about IMq. 

In order to study the first order quark-hadron phase transition it is necessary to specify the equation of 
state of the matter in both the quark and hadron states. So, at this stage, let us at first discuss the EoS for a 
first order phase transition |32j which in the quark phase can generally be given in the form 



Pq = 3a q T 4 +V(T), 
(ir 2 /90)g q , with g q = 16 + (21/2)N F - 



where a q 
potential, which is expressed as 



Pq = 

14.25 



a q T 



V(T), 
51.25 and N F 



V(T) = B + lT T 2 - a T T 4 



(23) 

2. V(T) is a self-interacting 

(24) 



where B is the bag pressure constant and, according to the results obtained in low energy hadron spectroscopy, 
heavy ion collisions and phenomenological fits of light hadron properties, the constants are given as B 1 ! 4 S 
(100 — 200)^^6^ [40], aT = 77r 2 /20 and 7t = m 2 /A with m s being the mass of the strange quark in the range 
m s e (60 — 200) MeV . The shape of the potential V reflects the physical model in which the quark fields are 
interacting with a chiral field formed out of the it meson and a scalar field. For the hadron phase, we choose a 
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cosmological fluid consisting of an ideal gas of massless pions and nucleons described by the Maxwell-Boltzmann 
statistics, having an energy density ph and pressure p/,, respectively. The equation of state is given by 

Ph{t) = \p h {T)^a^T\ (25) 

where — (ir 2 /90)gh and = 17.25. The critical temperature T c is defined by the condition p q (T c ) — Ph(T c ) 
[H], and is given by 



T = It + Vi't + 4(a g + a T - a*)B 
c V 2(a q + a T -a 7T ) ' 1 ' 

For m s = 20QMeV and B 1 / 4 = 200MeV the transition temperature is of the order of T c w 125MeV. A 
first order phase transition incorporates discontinuities for all the physical quantities like the energy density, 
pressure and entropy across the critical curve. The quark and hadron energy densities ratio at the critical 
temperature, p q (T c )/ p h (T c ), is of the order 3.62 for m s = 200MeV and B 1 / 4 = 200MeV. 



5 Dynamics of the universe during the first order transition 



After determining the equation of state, the physical quantities like the energy density p, temperature T, 
scale factor a(t) and Habble parameter H can be monitored during the phase transition. These quantities 
are determined by the Eqs. (|20|). ([22j) . ([23)) and ([25j) . Naturally, there are three stages to distinguish; before 
the phase transition where T > T c and the universe is likely to be in the quark phase, during the phase 
transition where T = T c and all the quark matter is changed to hadrons and finally, after the phase transition 
where T < T c and the universe is in the hadron phase. Now, using the equations of state of the quark 
matter, equation (|22j) . the variation of cosmological quantities such as the Hubble parameter, scale factor and 
temperature can be written as 



H(T) 



where after integration we have 



3a g — ax T 1 7t T 



3a„ 



T 



6o ? T 3 ' 



a(T) 



_ fl y(a r -3a ) 



3a„ 



exp 



1 It 1 
12 ~a q T i 



(27) 



(28) 



where ao is a constant of the integration and the variation of temperature for T > T c can be written as 



dT 
Ik 



rp3 



V 3a Q * 6a Q 



[(3a, - a T )T 4 + 7T T 2 + B + 2a\ 



6(3A- 1) 



(29) 



Figures [TJ [5] and |3] show the behavior of temperature as a function of the cosmic time r , Hubble parameter 
and scale factor as a function of temperature in a HL world filled with quark matter for different values of A 
For T = T c , the temperature and pressure are constant. In this case quantities like the entropy S = sa 
and enthalpy W — (p + p)a 3 are conserved. For later convenience, we replace p(t) by h(t), so that the volume 
fraction of matter in the hadron phase is given by 



p(t) = p h (t) + p q (l - h(t)) = p Q (l + mh(t)), 



(30) 



— 1. The beginning of the phase transition is characterized by h(t c ) 



where m 

scale and p(t c ) — pq, while the end of the transition is characterized by h(th 
signaling the end and corresponding to p{th) = Ph- To find the variation of hadron fraction we substitute 
equation (|30[) in the field equations and obtain 



= where t c is the time 
1 with th being the time 



dh 
~dt 



3(l + rh(t)) 



k 2 PQ 
6(3A-1) 



(1 +mh(t)) 



3(3A-1) 



1/2 



(31) 



where r = [pjj — Pq)/(pq + p c )- Figure @] shows variation of the hadron fraction h(r) as a function of r 
for different values of A. One may present variation of the Hubble parameter in this case as a function of 
the hadron fraction by using the field equation and equation (1301) , for which the result is shown in figure [S] 
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Figure 1: Variation of temperature in the T > T c region as a function of r = nt in a first order phase transition for different 
values of A : A = 3 (dotted-dashed curve), A = 2 (dotted curve), A = 1 (solid curve), A = 0.7 (long dashed curve), and A = 0.5 
(short dashed curve), we set -2j = 10 6 . 
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Figure 2: Variation of the Habble parameter in the T > T c region as a function of temperature in a first order phase transition 
for different values of A : A = 3 (dotted-dashed curve), A = 2 (dotted curve), A = 1 (solid curve), A = 0.7 (long dashed curve), and 
A = 0.5 (short dashed curve), we set -Sj- = 10 6 . 
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Figure 3: Variation of the scale factor in the T > T c region as a function of temperature in a first order phase transition for 
different values of A : A = 3 (dotted-dashed curve), A = 2 (dotted curve), A = f (solid curve), A = 0.7 (long dashed curve), and 
A = 0.5 (short dashed curve), we set = 10 6 . 
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Figure 4: Variation of the hadron fraction for T = T c as a function of r = nt in a first order phase transition for different values 
of A : A = 3 (dottcd-dashed curve), A = 2 (dotted curve), A = 1 (solid curve), A = 0.7 (long dashed curve), and A = 0.5 (short 
dashed curve), we set -% = 10 6 . 
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Figure 5: Variation of the Habble parameter in T 
different values of A : A = 3 (dotted-dashed curve), 
A = 0.5 (short dashed curve), we set = 10 6 . 



= T c case as a function of hadron fraction in a first order phase transition for 
A = 2 (dotted curve), A = 1 (solid curve), A = 0.7 (long dashed curve), and 
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Figure 6: Variation of temperature for T < T c as a function of r = nt in a first order phase transition for different values of 
A : A = 3 (dotted-dashed curve), A = 2 (dotted curve), A = 1 (solid curve), A = 0.7 (long dashed curve), and A = 0.5 (short dashed 
curve), we set = 10 6 . 
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Figure 7: Variation of the Habble parameter for T < T c as a function of temperature in a first order phase transition for different 
values of A : A = 3 (dotted-dashed curve), A = 2 (dotted curve), A = 1 (solid curve), A = 0.7 (long dashed curve), and A = 0.5 
(short dashed curve), we set = 10 6 . 

for different values of A. For t > th the universe enters into the hadronic phase. For the phase transition 
temperature of T c — YlhMeV we have pq = 5 x 10 9 MeV 4 and pn = 1-38 x 10 9 MeV 4 , respectively. For the 
same value of the temperature the value of the pressure of the cosmological fluid during the phase transition 
is p c = 4.6 x 10 8 MeV . Following the same procedure as that in the quark phase, that is, substituting 
equation (j2"5)l into (j2"2")l . the evolution of the scale factor can be obtained where a(rh) is now determined by 
using the conservation equation (|22[) and equations of state for T < T c which leads to a(T) = a(rh)T c /T. In 
the last stage of the first order phase transition, variation of temperature in the hadron phase (T < T c ) for a 
flat universe leads to 



(32) 



The variation of temperature in a first order phase transition as a function of r = kt and also variation of the 
Hubble parameter as a function of temperature are represented for different values of the coupling constant A, 
in figures |6]and[7l respectively. This parameter shows that small values of the coupling constant A reduce the 
temperature in these regions. Also, from figures [1] and [5] it is seen that for A > 1, the temperature is higher 
than in general relativity (A = 1) whereas for 1/3 < A < 1 the temperature is much lower. In the next section, 
we shall consider the evolution of temperature for a crossover phase transitions in Horava-Lifshitz theory. 




6 Crossover quark-hadron phase transition 

Indeed, our work is based on results obtained from two different sets of lattice calculations. The hot QCD 
collaboration analysis is based on two different improved staggered fermion actions, asqtad and p4, with a 
physical strange quark mass and physical u and d quark masses, m s /m Ut d = 10. It shows that the transition 
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region lies in the range T = (185 — 195)MeV, for recent published results see [37]. Results presented by 
the other collaborations Points to the value of the transition temperature in the range 150 — 170MeV, and 
changes with the observable used to define it [23J 23] . The most recent results for several physical quantities 
are presented in [12] with the equation of state at N t = 6, 8, 10 and 12 given in [15]. We shall rely on the above 
two cases of analysis to determine the EoS in the context of the crossover phase transition in this work. In the 
first case, the equation of states in both phases are given by the hadron resonance gas (HRG) model at low 
temperatures and lattice QCD data at high temperatures, which are combined with a simple parametrization 
given in [26]. The equation of state has been calculated with the so-called p4 and asqtad actions on lattices 
with temporal extent N T = 4, 6 and 8 [Ml 1551 157] . 

Calculation of the thermodynamical observables in lattice QCD proceeds through the calculation of the trace 
of the energy momentum tensor 0(T) — e — 3p, also known as the trace anomaly. Different thermodynamical 
observables can be obtained from the trace anomaly through integration. The pressure can be written as 

p(T) p(T ) f T dT' 

;(p-3p). (33) 



T 4 Tq Jrp T' 5 

Here To is an arbitrary temperature value that is usually chosen in the low-temperature regime where the 
pressure is exponentially small. This method is known as the integral method [38] . Finite temperature lattice 
calculations are mostly performed at a fixed temporal extent N T with the temperature being changed by 
varying the lattice spacing a, (T = l/N T a). Therefore, the lattice spacing becomes smaller as the temperature 
is increased. As a consequence, the trace anomaly can be calculated in the high temperature region accurately, 
whereas at low temperature is affected by the large discretization effects. Thus, to construct a more realistic 
equation of state one should use the lattice data for the trace anomaly in the high temperature region, T > 
250MeV, and the HRG model otherwise, T < ISOMeV. 

The HRG model with modified masses is believed to describe the lattice data rather well up to temperatures 
of about l80MeV. In the middle region, 180MeV < T < 250MeV, the HRG model cannot be trusted, whereas 
discretization effects in lattice calculations could be large. Therefore, at high temperatures the trace anomaly 
can be parameterized in a polynomial form in general, see e.g. |37) . We shall thus use the following ansatz for 
the high temperature region 

~ 3g) = ^2_,^,_£l_,_£2_ ( oa) 

This form is shown to work well in the temperature range of interest and is flexible enough to be matched 
to the HRG result in the low temperature region. The resulting parametrization of the trace anomaly and 
QCD equation of state obtained using these requirements are labeled as s95p — vl, s95n — vl and s90/ — vl, 
as shown in table 1. The labels "s95" and "s90" refer to the fraction of the ideal entropy density reached at 
T = 800MeV (95% and 90%, respectively). The labels p, n and / refer to a specific treatment of the peak of 
the trace anomaly or its matching to the HRG. Values of the parameters To, d 2 , ^4, ci, C2, n\ and n 2 in are 
given in table 1. The following form shows a parameterized form for the HRG 

= aiT + a 2 T 3 + a 3 T 4 + a 4 T 10 , (35) 

where a x = 4.654GeV- 1 , a 2 = -879GeV- 3 , a 3 = 8081G e y- 4 and a 4 = -7039000G e y- 10 . The importance of 
the above parametrization and integral method is that we can obtain the thermodynamical quantities for the 
low temperature region as 

p(T) = (a + /3 )T 4 + ai T 5 + ^a 2 T 7 + ^a 3 T 8 + ^T 14 , (36) 

p(T) = 3(a + /3 )T 4 + 4 ai T 5 + 2a 2 T 7 + ? -a 3 T 8 + ^a 4 T 14 , (37) 

where o?o and /3q are the value of integration in the lower limit and the value of Pq/Tq when To = 70MeV, 
respectively. For the high temperature region, the pressure and energy density are given by 

p(T) - -~di - —T 2 + (a, + A)T 4 - flT 4 "™ 1 - ^ T 4 ' n \ (38) 
4 2 n\ n 2 

p{T) = 3(a, + ft)T 4 - \d 2 T 2 + ]d 4 + ^(ni - 3)T 4 -' 11 + ^(n 2 - 3)T 4 "" 2 , (39) 
2 4 m n 2 
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Table 1: The values of parameters for different fits to the trace anomaly I26| . 





d 2 (GeV) 


d^GeV^) 


Cl {GeV ni ) 


c 2 (GeV n2 ) 


m 


n 2 


To(MeV) 


s95p 


0.2660 


2.403 x 10" 3 


-2.809 x 10" 7 


6.073 x 10- 


23 


10 


30 


183.8 


s95n 


0.2654 


6.563 x 10~ 3 


-4.370 x 10~ 5 


5.774 x 10" 


-6 


8 


9 


171.8 


s90f 


0.2495 


1.355 x 10" 2 


-3.237 x 10" 3 


1.439 x 10- 


14 


5 


18 


170.0 



Table 2 : The values of parameters appearing in equation l|40|l describing the trace anomaly for (jf =2 + 1 and ri^ = 2 + l + l 

JUL 





ho 


hi 


h 2 


fo 


h 


h 


9i 


92 


2 + 1 


0.1396 


-0.1800 


0.0350 


2.76 


6.79 


-5.29 


-0.47 


1.04 


2 + 1 + 1 


0.1396 


-0.1800 


0.0350 


5.59 


7.34 


-5.60 


1.42 


0.50 



where i = 1,2,3 corresponds to different parameterizations which are given in table 1. So, on and ft refer to 
the values of integration in the lower limit and the value of Po/Tq in each case. Equations of state are also 
available in a tabulated form in [36] . 

As was mentioned above, there is another set of data which we will use for determining the EoS in the 
present work for which a simple parametrization of the trace anomaly, given for rif = 2 + 1 flavor and also for 
an estimate of rif = 2 + 1 + 1 flavor equation of state, can be found in [45]. We notice that these results are 
obtained for Symanzik improved gauge and a stout-link improved staggered fermion action on lattices with 
temporal extent N T = 6,8 and 10. The thermodynamical quantities are presented as functions of temperature 
in the range 100 • • • lOOOMeT^. A global parametrization of the trace anomaly as a function of temperature is 
given by the following fitted function 



T 4 



exp (—hi/t — h 2 /t 2 



h 



/o(tanh(/ 1 t+/ 2 ) + l) 
l+git+g 2 f 2 



(40) 



where the dimensionless t variable is defined as t — T/(200MeV). The parameters can be found in 2able 2. In 
this case, the numerical solution is used for determining the EoS. The results of this calculation are presented 
in the following section. 

7 Dynamical evolution in the crossover phase transition 

In this section we shall study the physical quantities of interest, that is, energy density p and temperature 
T. For the first set of lattice data [37] mentioned before, these quantities are obtained from the Friedmann 
equation ([20]). conservation equation ([22]) and the equations of state determined by equations ([36]) and (]37[) 
for low temperature region and equations (]38j) and (|39[) for high temperature region in the crossover phase 
transition. By assuming a first order phase transition, these quantities can be obtained from the Friedmann 
equation ([20]) . conservation equation ([22]) and equations of state ([23]) and (1251) . determined for T > T c and 
T < T c , respectively. Using the conservation equation d22l and equation of state which is used for describing 
the high temperature region, we have 



a 
a 



(41) 



12(q, + ft)T 3 - d 2 T + ^(m - 3)(4 - n^T 3 "" 1 + %(rc 2 - 3)(4 - n 2 )T 3 ' 
3(-d 2 T 2 + 4(a, + ft)T 4 + £(m - 4)r 4 -™! + £(n 2 - 4)7^) 



Substituting equation ([41]) in equation ([20]) , the basic equation describing the evolution of temperature of the 
universe in the high temperature region can be written as 

dT _ 3(-d 2 T 2 + 4( ai + A)T 4 + %( ni - 4)T 4 -"i + ^(n 2 - 4)T 4 -"^) 

dt ~ 12(a, + &)T 3 - d 2 T + £(ni - 3)(4 - n^T^ + ^(n 2 - 3)(4 - n 2 )T^ 



, 3(a, + ft )T 4 - \d 2 T 2 + \d A + £L ( ni _ 3 )T 4 -"i + a („ 2 _ 3 )T 



4-n 2 



6(3A- 1) 



+ 



kA k 2 B k 3 C 



3(3A-1) a 2 a 4 (i) a 6 (i) 



(42) 
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Figure 8: Variation of temperature in the high temperature region as a function of r = nt for the first parametrization 
corresponding to To = !70.0MeV . The behavior of temperature is shown for different values of A : A = 3 (dotted-dashed curve), 
A = 2 (dotted curve), A = 1 (solid curve), A = 0.7 (long dashed curve) and A = 0.5 (short dashed curve), we set = 10 10 . 



where 



.4 



2" 



(3A-1) 



D 



4(37! +72) c _ 8(9cti + 3<j 2 + <t 3 ) 



(3A - 1) 



(3A-1) 



(43) 



and other parameters are determined as in the previous section. The coefficients of the curvature are absent 
in a spatially flat universe (k = 0), thus the coefficients A, B and C are zero, so that equation (14"2"T) reduces to 

dT _ 3(-d 2 T 2 + 4(a, + &)T 4 + %-(m - 4)T 4 -»i + %(n 2 - 4)r 4 -"^) 

dt ~ 12(a l + /3 l )r 3 -d 2 T+^(n 1 -3)(4-n 1 )T3""i + ^(n 2 -3)(4-n 2 )T3-^ 

V(3(q a + &)T 4 - \d 2 T 2 + |d4 + f^i - i)T^ + %(n 2 - 3)T 4 -^) + 2a 

6(3A- 1) 



(44) 



which has been modified relative to the general relativity result with parameter A. 

Moving forward, we first study the first order phase transition where the curvature term is zero. We 
then focus attention on the study of a crossover phase transition with the same conditions. The variation of 
temperature as a function of r = ret in the high temperature region is presented for a flat universe for different 
values of the coupling constant A, in figure [H The behavior of this parameter shows that at small values of 
the coupling constant A, the temperature is reduced in this era. As can be seen in these figures, for A > 1 the 
temperature is higher than the general relativity case where A = 1, whereas for 1/3 < A < 1 the temperature 
is much lower. 

For the low temperature region, T < 170MeV, using the conservation equation (f22j) and equations of state 
(l38l) and (l39|) and also equation (l20|) , the basic equation describing the evolution of temperature of the universe 
can be written as 



dT 
~dt 



12(q + Pq)T + 15qiT 2 + 7a 2 r 4 + 6a 3 r 5 + fa^T 11 



12(a " 
2 3(a 



Ai) + 20aiT + 14a 2 T 3 + 14a 3 T 4 
- (3 Q )T 4 + 4 &1 T 5 + 2a 2 T 7 



i«3 



^8 



CL4T 



14 



6(3A- 1) 



(3A-1) 



(45) 



where a dot represents derivative respect to time. The variation of temperature as a function of parameter 
r = nt in the low temperature region is represented, in the case of fc = 0, for different values of the coupling 
constant A, in figure [9] The behavior of this parameter shows that small values of the coupling constant A 
reduce the effective temperature. Figure [9] shows that for A > 1, the temperature is higher than in general 
relativity (A = 1) whereas for 1/3 < A < 1 the temperature is much lower, again relative to general relativity. 

At this stage, we turn to the second set of lattice data [45] mentioned before and repeat the above calcu- 
lations. However, this time there is no analytical expression for the equation of state and one has to resort 
to numerical calculations to obtain the evolution of temperature for the parameterization given in equatiori40l 
Variation of temperature for rif = 2 + 1 and rif = 2 + 1 + 1 flavors for different values of A is presented in Figs. 
1101 and 1111 respectively. One realizes that the effect of A on the variation of temperature in the early universe 
in this case is the same as the results obtained by the "hot QCD" collaboration. 
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Figure 9: Variation of temperature of the low temperature region as a function of r = nt for different values 
of A : A = 3 (dotted-dashed curve), A = 2 (dotted curve), A = 1 (solid curve), A = 0.7 (long dashed curve), and 
A = 0.5 (short dashed curve), we set = 10 6 . 
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Figure 10: Variation of temperature for nj = 2 + 1 flavor as a function of r = nt in a smooth cross-over phase transition for 
different values of A : A = 3 (dotted-dashed curve), A = 2 (dotted curve), A = 1 (solid curve), A = 0.7 (long dashed curve), and 
A = 0.5 (short dashed curve) where we set = 10 6 . 



200 - 




100 u — 1 — 1 1 — 1 — 1 1 — 1 1 — 1 — 1 1 — 1 1 — 1 — 1 1 — 1 — 1 1 — u 

2.X10" 6 4.x 10- 6 6.X 10- 6 8.x 10- 6 0.00001 

t MeV -2 



Figure 1 1 : Variation of temperature for rif = 2 + 1 + 1 flavor as a function of r = nt in a smooth cross-over phase transition for 
different values of A : A = 3 (dotted-dashed curve), A = 2 (dotted curve), A = 1 (solid curve), A = 0.7 (long dashed curve), and 
A = 0.5 (short dashed curve) where we set = 10 e . 
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8 Conclusions and Remarks 



In this work, we have studied the QCD phase transition in the context of Horava-Lifshitz gravity in the early 
universe, where differences from what is predicted by general relativity theory are likely to be noticeable. We 
investigated the effects of the coupling constant on the evolution of physical quantities, relevant to the physical 
description of the early universe in the case of first order and smooth crossover phase transitions. We used 
phenomenological equation of state which led to first order phase transition and investigated the influence 
of the coupling constant on the variation of temperature and Hubble parameter. The results show that the 
Hubble parameter decreases with time and that an increasing coupling constant decrease the value of this 
parameter before and after the phase transition. In addition, higher values of the coupling constant increase 
the effective temperature of the quark and hadron phases and also decrease the hadron fraction in mixed phase. 
In this work, two different sets of lattice simulation data which led to two different cases for the equation of 
state with smooth crossover phase transition, are used. In the first case, it led to the energy density and 
temperature in both high (T > 250MeV) and low (T < 170MeV) temperature regions with the equation of 
state obtained by the "hot QCD" collaboration which is believed to be responsible for the emergence of hadrons 
from the quarks and gluons. In the second case, we used results obtained with the lattice data provided by the 
Wuppertal-Budabest collaboration. Comparison of these two different sets of results shows good agreement 
for the influence of the coupling constant on the variation of temperature in the context of the HL gravity. 
We found that the temperature evolution of the early universe in Horava gravity is different from the idealized 
standard FRW model, so that for a small value of the coupling constant A, phase transition occurs and results 
in describing the effective temperature in both high and low temperature regions. We also noted that the 
general behavior of the physical quantities mentioned above would not differ markedly if a first order phase 
transition is assumed to have occurred. 
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